A geometric approach to the diophantine Frobenius problem 
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ABSTRACT. It turns out that all instances of the diophantine Frobenius problem for three coprime aj 
have a common geometric structure which is independent of arithmetic coincidences among the a^. By 
exploiting this structure we easily obtain Johnson's formula for the largest non-representable z, as well 
as a formula for the number of such z. A procedure is described which computes these quantities in 
O(log(maxai)) steps. 



1. INTRODUCTION 

For an n-tuple a = (oi, 02, . . . , a„) of positive integers we denote by T := T{a.) the set 
of natural numbers z that can be written in the form 



^XfcOfc, Xfc G N := {0,1,2...} , 



fc=i 



and by F := -F(a) the set of natural numbers z that cannot be so represented. If 
gcd(ai, . . . ,an) = 1 then it is easily seen that all sufficiently large numbers z are in 
T. It follows that in this case F is a finite set, so there is a largest non-representable 
^ ' number g{a) := maxF. To compute this number and maybe even the cardinality A'^(a) 

Cp , of -F in terms of oi, . . ., an constitutes the so-called diophantine Frobenius problem. 

We recommend [4], printed in 2005, as a comprehensive source of material about this 
problem; the bibliography alone contains about 500 items. 

^^ I The case n = 2 was first considered and solved by Sylvester [6], [7]. He proved: 



Proposition 1. If oi, 02 are > 1 and coprime then 



. J ■ (ai — l)(a2 — 1) 

P>^ ■ 5(01,02) = 01^2 -«! -a2 , N{ai,a2) = • (1) 



The proof of this result follows from inspection of Fig. 1 and is given at the beginning 
of the next section. 

This paper deals with the case re = 3. We shall give a natural geometric description of 
the set F from which an explicit answer to the Frobenius-3-problem can immediately be 
read off. In order to formulate our result we introduce the quantities 

I3 :=min{/ G N>oUa3 eT(ai,a2)} O . 

Here and in the sequel the sign O indicates that there are three such formulae in all, 
whereby the other two are obtained by cyclic permutation 1— )-2— )-3— )-lof the indices. 
About the li the following can be said right away (cf. [2], Theorem 3 and eq. 26): 



Proposition 2. Assume that the three numbers oi, 02, 03 are pairwise prime and that 
^i^2 (l<i<3). Then the minimal representation 

ka3= X3iai+X32a2 , X31,X32GN O (2) 

of liGi (1 < i < 3) is uniquely determined, and one has 

Xij > 1 (all i^j) . (3) 

Furthermore the li are coupled to the Xij through 

k = a;i3 + X23 O . (4) 

We now state our main result; it will be proven in section 3: 

Theorem 3. Assume that the three numbers oi, 02, 03 are pairwise prime and that 
li>2 (1 < i < 3). Then 

g{a.) = li I2 h + max{xi2X23a;3i, a;2ia;32a;i3} - V!. Oi ; (5) 

A^(a) = Vj](/i-l)a, - /1/2/3 + lV (6) 

The assumption li > 2 means that none of the a^ is "superfluous". If, e.g., ^3 = 1 then 
F{a-i,a2,as) = F{ai,a2); this case is handled in Proposition 1. When the given a^ are 
not pairwise prime then there is a way to get rid of common factors, see [2], Theorem 2. 
Our formula (5), resp. its preliminary version (8), appears as Theorem 4 in [2] and on 
p. 35 of [4]. Note, however, that the proof given in [4] uses heavy algebraic machinery 
and is deferred to a later chapter. 

Example 1. Let ai := 2n — 1, 02 := 2n, 03 := 2n + 1 for an n > 2. As 02 = 1 and 03 = 2 
(mod ai) the smallest multiple of ai in T(a2, as) is 02 + (n — 1)03; similarly the smallest 
multiple of 03 in T(ai,a2) is noi + 02, and obviously the minimal representation of 02 
is 2a2 = 01+03. Altogether we have 

h =n + l, X12 = 1, a;i3 = n - 1; I2 = 2, a;2i = 1, 0:23 = 1; ^3 = n, a;3i = n, ^32 = 1 ; 

so Theorem 3 gives 

g{a) = 2n{n + 1) + max{n, n — 1} — 6n = 2n — 3n , 

N{a) = - (n(2n - 1) + 2n + (n - l)(2n + 1) - 2n(n + 1) + l) = n^ - n . 

For the computation of the li and the Xij we propose the so called Lagrange algorithm 

- a kind of two-dimensional euclidean algorithm modeled after a Gram-Schmidt-process 

- which takes 0(log(maxai)) steps. The resulting procedure is developped in sections 
4 and 5 of this paper. In [4] several other algorithms for ^(a) are described, notably the 
algorithm of R0dseth [5] which is an improved version of an earlier continued fraction 
algorithm by Selmer &: Beyer. 
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Fig. 1 



2. PRELIMINARIES 



Proof of Proposition 1. We draw in the (xi,X2)-plane the directed graph F with vertex 
set Z^ and edges of unit length connecting neighboring lattice points in the direction of 
increasing xi resp. X2, see Fig. 1. For given oi, 02 G N>i we define the height function 

f{xi,X2) := flixi +02X2 . 

Two points in Z^ have the same height iff they differ by a vector u S L := Z^ n /~^(0). 
Since ai and 02 are coprime the set L is the one-dimensional lattice formed by the 
vectors u^ := {ka2, —kai), k £ Z. 

An integer z G N can be represented in the form z = xiai + X2 02 with xi, X2 G N iff 
there is a directed edge path in F connecting a point uj. £ L with a point x G Z^ of 
height /(x) = z. Now the lattice points that can be reached from a given u^ G L lie in 
the set Qk '■= u^ -(- M>q , a first quadrant with origin at u^, and the set of all possible 
end points of such paths consists of the lattice points in the union Q := |J^ Qk of these 
quadrants. 

The lattice points of positive height that cannot be reached from one of the points u^ 
are the interior lattice points of the rectangular triangles A^ with vertices u^, Ufc_|_i, 
Ufc + (a2, 0). The largest occurring height in such a A^ is given by the first formula (1), 
and using symmetry one concludes that each A^ contains exactly (oi — l)(a2 — l)/2 
lattice points in its interior, which all have different heights. D 

For later purposes we note the following: Any lattice point in Aq can be connected by 
an admissible path to the point (02, 0) of height 0102- It follows that a number z > is 
in F{ai,a2) iff there are integers ki, k2 > 1 such that z = aia2 — kiai — ^202- 

Proof of Proposition 2. First we show (3). Assume, e.g., that X13 = 0. Then we have 
liai = Xi2a2, and as oi, 02 are coprime it follows that li > 02. On the other hand, from 
^3 > 2, i.e., 03 G F(ai,a2) it follows that there are integers ki, k2 > 1 with 

03 = aia2 — kiai — k2a2 = (02 — ^i)ai — ^202 • 



Whence we would have (02 — ki)ai = ^202 + Q3 which would imply /i < 02 - a con- 
tradiction. 

We next show (4). Let the Xij be determined such that (2) holds. Then the quantities 

/U3 := ^3 - 2;i3 - 2^23 O 

satisfy jUiOi + jU2a2 + ^303 = 0. If the Hi do not all vanish then up to a permutation of 
the tti we must have one of the following: 

(a) ^1 > , ^2 < , ^3 < ; (b) ^1 > , ^2 > , ^ug < . 

In case (a) it follows that 

(/i - a;2i - X3i)ai = ^lOi = (-/i2)a2 + (-/U3)a3 , 

contradicting the definition of /i. In case (b), from 

(xi3 + X23 - ^3)03 = -^303 = MiOi + ^202 

it follows by definition of ^3 that X13 + X23 > 2/3, whence, e.g., X13 > ^3. By definition 
of the Xij we now have the representation 

ih - X3i)ai = (a;i2 + a;32)a2 + (a;i3 - ^3)03 

which again contradicts the definition of li. 

As (4) is true for all possible choices of X21 and X31 consistent with the definition of I2 
and ^3, and as these choices can be made independently, it follows that there is in fact 
no choice at all, which means that the Xij are indeed uniquely determined. D 

3. PROOF OF THE MAIN RESULT 

We now come to the proof of Theorem 3. Inspired by the proof of Sylvester's result for 
n = 2 we embed the problem into the following geometric setup: Consider the integer 
lattice Z^ in euclidean (xi, X2, X3)-space M^. We use Z^ as vertex set of a directed graph 
r whose edges are the segments of unit length connecting neighboring lattice points in 
the direction of increasing xi, resp. X2, X3. Given oi, 02, 03, we again define the height 
function 

/(xi, X2, X3) := ai xi + 02X2 + 03 X3 

which on the one hand is just a linear functional on M^ and on the other hand assigns a 
height /(x) to each lattice point x G Z^. The kernel H := /~^(0) of / is a plane through 
the origin of R^ and contains the Frobenius lattice L := H nZ^ of integer solutions to 
the equation /(x) = 0. 

Lemma 4. (a) Let niiai +771202 = 1 with rrii G Z. Then the vectors ei := (02, — ai,0), 
^2 '■= (o3"^i) 03"J2; — 1) form a basis of L. 

(b) The fundamental domain of the lattice L, when projected to the plane Xi = 0, has 
area ai (1 < i < 3). 



Proof, (a) One easily checks that /(ei) = /(e2) = 0, which means that ei, e2 € L. 
On the other hand, let u = {ui,U2,U3) be an arbitrary point of L. Then u + 1^362 = 
{u'i,u'2,0) G L which implies aiu'i + 021*2 — 0- Since oi, 02 are coprime it follows that 
{ui,u'2,0) = kei for a /c € Z, whence we have u = kei — ^362. 

(b) It suffices to compute the vector product 



ei X e2 = (ai,a2, (miai + 771202)03) = (01,02,03) . 



n 



The three lattice vectors 

fi := (^1,-2:12,-2:13), h ■■= (-2:21,^2,-2:23), f3 := (-2:31,-2:32,^3) e L 

encoding the data li, Xij will play a special role. We shall call any vector of the form f^ 
or — fi a basic vector and any set of three essentially different vectors among the =b f^ a 
solution basis for the Frobenius problem at hand. 




Fig. 2 



The connection of the three-dimensional structure described so far with the diophantine 
Frobenius problem is the following: A natural number z is in T'(a) iff there is a lattice 
point u G L and a directed edge path in F beginning at u and ending in a point x G Z'^ 
of height /(x) = z. Now the lattice points that can be reached from a given u G L lie 
in the set Ou := u + M>q , an octant with origin at u, and the set of all possible end 
points of such paths consists of the lattice points in the union Q := IJ^gL ^u °f these 
octants. It follows that r(a) = f{1,^ n Q). The boundary dO, is L-periodic; it consists 
of three L-shapes per octant and touches the plane H in the points of L. It looks like 
a washboard and is depicted in Fig. 2, as seen from below. The proof of Theorem 3 
essentially consists in understanding this figure. 



Consider, e.g., the positive xs-axis. It is an edge of the octant Oq and belongs to the 
boundary of 17 until it is intercepted at p = (0, 0, us) by a wall xs = const, belonging to 
another octant Ou, u = {ui,U2,U3) G L. For this to happen it is necessary that ui < 0, 
U2 < 0, Us > 0, and that there is no point u' = (it'^, U2, ^3) G L with u'l < 0, u'2 < and 
U3 < Us- This means 

U3 = minju e N>o I 3ui,U2 € Ijko: oiui + 02^2 + 03^ = 0| ^ ^ 

"1 (7) 

= minju e N>o | naa G T{ai,a2)} , 

from which we deduce 

U3 = h , -Ul = X31 , -U2 = X32 , 

i.e., {ui,U2,U3) = {3 ; and similarly for the other li, Xij. 

The lattice points x of positive height that cannot be reached from a point u G L are 
the interior lattice points contained in the region W enclosed between H and dQ. As 
seen in the figure, the restriction / \ dQ takes local maxima at the points qi, q2, and 
the interior lattice points of maximal height are qi — (1, 1, 1) or q2 — (1, 1, 1) and their 
equivalents modL. It follows that the maximal non-realizable height g{a) is given by 

g{a.) = max{/(qi),/(q2)} - ^ a^ = ^3 03 + max{x2iai,Xi2a2} -'^."-i ■ (8) 

Since the three L-shapes have areas a^ by Lemma 4(b), we deduce from Fig. 2 that the 
Gi satisfy 

ai = a;i2^3 + a;i3a:;32 , 02 = a;2iZ3 + a;23a;3i , 03 = ^1^2 - a;i2X2i O . (9) 

Substituting these expressions into (8) one arrives at the symmetric formula (5): 

5(a) = I1I2I3 -^32:123:^21 +max{x2iXi2/3 + a:2ia:;i32;32,xi2X2iZ3 + a:i2a;232;3i} -y^^.o-i 

= ^i^2^3 +max{x2lXi3X32,Xi2X23X3i} - V" 



'I 



We now come to the proof of formula (6). We have to count the number zq of interior 
lattice points in the quotient W := W/L. As W does not have a simple description 
in terms of inequalities we are going to determine zq "from the outside" by means of 
a three-dimensional analog of Pick's area formula. Let zi denote the total number of 
relative interior lattice points in the three L-shapes; similarly, let Z2 be the total number 
of relative interior lattice points on the reentrant edges of W and Z3 be the number of 
such points on the protruding edges of W. Then we have the following formula: 

13 17 

Lemma 5. vol(VF) = zq + -zi + -Z2 -|- -23 + t , 

whereby the last term incorporates the contribution of the corners of W . 

Proof. We perform a "Gedankenexperiment" used already in [1] for a proof of Pick's 
area formula. Imagine that at time a unit of heat is concentrated at each point of Z'^. 



This heat wih be distributed all over space by heat conduction, and at time oo it will be 
equally distributed in space with density 1. In particular, the amount of heat contained 
in W will be vol(VF). Where does this amount of heat come from? For symmetry reasons 
there is absolutely no flux across the unit squares of d^, and, again by symmetry, the 
net flux across H/L is as well. As a consequence, the final amount of heat within W 
comes from the interior lattice points, counted by zq, and from the lattice points on the 
boundary of W. The lattice points counted by zi send half their heat into W, whereas 
the corresponding factor is | for the points counted by Z2 and j for the points counted 
by Z3. Furthermore W possesses two protruding corners qi and q2 which contribute 
I each, three "L-corners" contributing |, and finally the reentrant corner on H which 
contributes | as well. D 
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Fig. 3 



In order to compute vol(W^) directly we use the L-shape A in the plane 3:3 = as 
fundamental domain. Fig. 3 shows A, as seen from the positive xa-direction. The 
vertical prism K determined by A and the plane 



H : 



X3 



1 
03 



(oixi + 022:2) 



is a representative for W. To compute the volume of K we split A into two rectangles 
and work with the heights of K in their midpoints. We obtain 



uim 2:31X12/ X31 2x32 +X12 

vol{W) = Ui — - + a2 

03 V 2 2 



hx32 f h , 3:32 \ 
^^(ai- + a2— j 



which using (2) and (9) can be brought into the following symmetric form: 

3l(^) = ^ (I]^ ^i «. -^1/2/3' 



vol 



(10) 



We now compute the quantities zi, Z2 and z^. - The horizontal L-shape A in Fig. 3 
has area a^ by Lemma 4(b) and 2(?i + I2) lattice points on its boundary. Therefore by 
Pick's area formula for the plane the number of interior lattice points on A is given by 
fla ~ (^1 + h) + li and we obtain 
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:iii 



Inspection of Fig. 2 shows that 

Introducing (10), (11) and (12) into Lemma 5 we get 

uw^ 13 17 

zo = vol[W) - -zi - -Z2 - 4^3 - I 



(12) 



^ (5]x^^ - i)«i " ^1^2/3 + 1) ■ 



This concludes the proof of Theorem 3. 



D 



4. FINDING A SOLUTION BASIS 

In order to make Theorem 3 useful we have to establish a procedure to compute the 
quantities li, Xij. Our argument takes place in the plane H. Fig. 4 shows H as seen from 
the tip of the vector a, the points of L are again marked by bullets. The three planes 
Xi = intersect H in three lines gi through the origin which altogether divide H into six 
sectors of various widths. The line (73, spanned by the vector V3 := (— 02, ai, 0), is at the 
same time a level line of the linear function X3 restricted to H. Equation (7) can now 
be interpreted as follows: In order to find Is we have to translate the line §3: X3 = in 
the direction of increasing 2:3 (marked by an arrow in Fig. 4) until it hits for the first 
time a lattice point in the sector a;i<0 A X2 <0. The lattice point obtained in this way 
is the point £3. Translating similarly the lines gi and 52 one obtains the lattice points 
fi, £2 in the appropriate sectors. 



9l- 2:1 = 




92- ^2 = 



^3' ^3 



a;Q = 



Fig. 4 



Being minimizers of some sort the fi £ L tend to be short. Now there exists a revered 
algorithm (attributed to Lagrange, Gauss and others, see [3]) which finds the shortest 
vector u of the lattice L, and we plan to make use of this algorithm. But if u happens 
to lie in a sector of width > 60°, as in Fig. 4, there is no guarantee that u coincides with 
the basic vector f^ (or — f^) in that sector. For this reason we change the metric in such 
a way that the three lines Xi = intersect at angles of 60°. 

Lemma 6. For a suitable scalar product (x, y) := x Qy the three lines Xi = in H 
intersect at angles of 60°. 

Proof. The three directions in question are 

Vi := (0,-03,02) , V2 := (03,0,-01) , V3 := (-02,01,0) . 

Consider now the linear map P : M'^ — )• M^ given by the matrix 



P :-- 



-oi 02 

V303 



The kernel of P is spanned by the vector (02, oi, 0) ^ H, therefore the restriction P \ H 
maps H bijectively onto the euclidean plane E := M?. One easily computes 



1 V3n 



Pvi = 20203! -^, ^) , Pv2 = 20103 (--, -^ ) , Pv3 = 20102(1,0), 



1 

2' 



2 



which shows that in the image plane the lines Qi intersect at angles of 60°. 

Pulling back the euclidean scalar product in i? to if one obtains there the new scalar 
product 

(x, y) := (Px)T Py = xTpT p y = xT Q y , 

where the integer matrix Q := P P is given by 



Q 



ai 


-01O2 





01O2 



oi 




3oi 



n 



In what follows, | • | denotes the norm corresponding to the scalar product (•,•). 

Lagrange's algorithm (to be described in the next section) produces a reduced basis 
(u, v) of the Frobenius lattice L. This means that u is a shortest nonzero vector in L 
and that v is a shortest vector in L \ Zu; in particular, |u| < |v|. 

Theorem 7. Under the hypotheses of Theorem 3, let (u, v) be a reduced basis of the 
Frobenius lattice L. Assume that ui < 0, U2 < 0, n3 > and put A := V3/U3. Then the 
three vectors 

u , v_ := V — [A] u , v_|_ := V — [AJ u 

form a solution basis for the Frobenius problem determined by the o^. 

9 




53 



Fig. 5 

Proof. We argue in the (x, y)-plane E, but omit the P in our notation. The hues gi 
enclose angles of 60°, creating (closed) sectors 5*^- (1 < A; < 6), see Fig. 5. After scaling 
we have u = (— sin^,cos^) with \6\ < ^, whence u G 5*2, and we may assume that the 
lattice line ^ || u containing v is to the right of u. We shall show that the points u, v_ 
and v+ are basic points in the sectors Sq, Si and 5*2 respectively. 

We begin with the remark that in fact |^| < vr/G. Assume to the contrary that, e.g., 
6 = 7r/6. Then u = (— 1/2,\/3/2) G gi, meaning ui = 0. In this case u could not be 
basic by (3). It follows that there would have to be a lattice point in ^2 with y-coordinate 
< \/3/2. But there is no room for such a point since the interior of the unit circle is 
forbidden. 

Let b := V — Au = (6,0) be the point where I intersects g^. Then b = h/ cos9 where h 
denotes the distance from to i, whence h > \/3/2. We write (cos (/>, sin 0) =: e^. 

Lemmas, (a) (e^/g, u) < (e^/6,b) , (b) (e.^/g, -u) < (e_^/6,b) . 

Proof. The left sides of (a) and (b) are sin(-| — 0) and sin(^ + 6*) respectively, so they 

both are < sin(^ + |^|). On the other hand the right sides of (a) and (b) both have the 

same value 

vr TT h 3 

cos — = cos > . 

6 6 cost* 4 cosy 

It remains to prove that for < 9 < ^ one has 



/TT 

2sin - 

Ve 



cos < - . 



But here the left side can be written as sin(-|- + 2^) + sin(-|-) which is < 1 + ^. D 

Put c_ := b — u, c+ := b + u. Then from Lemma 8(a) it follows that {e^/Q,C-) > 0, 
whence c_ G int(S'6), and analogously from Lemma 8(b) it follows that (e_^/6)C+) > 0, 
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whence c_|_ € int(S'i). From this we conclude that the lattice points v_ G [b,c_] and 
v+ e [b, c+] lie in Sq and 5*1 respectively. 

As c_|_ G int(S'i) the line i intersects 52 at a y-level > cos 9, whence all lattice points on 
£n ^2 have a larger y-level than u, and similarly, as — c_ € int(S'3), the line i' intersects 
gi at a y-level > cos 9, whence all lattice points on i' n S2 have a larger y-level than u. 
This implies that the vector u is basic in its sector 5*2. 

Note that {e„/Q,u) > 0, whence going upwards along if] Si the distance to yi increases. 
Since v_|_ is the first lattice point met along this path, v_|_ is basic for the sector 5*1, 
unless there were an even better lattice point on the parallel to i through the point 2b. 
But the latter is prohibited by the inequality (e^r/g, c_|_) < {e„/Q, 2b) which follows easily 
from Lemma 8(a). 

Similarly one has {e_^/Q,u) < 0, and this implies that going downwards along i f] Sq 
the distance to 52 increases. Since v_ is the first lattice point met along this path, v_ 
is basic for the sector Sq, unless there were an even better lattice point on the parallel 
to i through the point 2b. But the latter is prohibited by the inequality {e_^/Q,c_) < 
(e_7r/6)2b) which follows easily from Lemma 8(b). D 

5. LAGRANGE'S ALGORITHM 

Lagrange's algorithm, as it is called in [3], takes an arbitrary basis (u, v) of the Frobenius 
lattice L as input and in a certain number of steps arrives at a reduced basis of L. An 
essential accessory to the calculations is the Gram matrix 



G:=G(u,v) := 



(u,u) (u,v) 
(v,u) (v,v) 



of the current basis (u, v); it is updated along with the basis vectors. 

The following box is taken from [3]. The subscript < to a basis indicates that one 
assumes |u| < |v|, and [ • ] denotes the nearest integer function. 



Input: A basis (u, v)< with its Gram matrix G = {gij)i<i,j<2- 
Output: A reduced basis of L with its Gram matrix. 




1. Repeat 




2. r := V — xu, where X := 51,2/51,1 • 

When computing x, also compute the remainder y := yi,2 — a;yi,i. 

3. v:=u. 




4. u:=r. 




5. Update the Gram matrix as follows: swap 52,2 and yi,i; then let yi,2 
51,1 :=5i,i -x{y + gi^2)- 


= y and 


6. Until |u| > |v| . 




7. Return (v, u)< and its Gram matrix (setting 52,1 := 51,2)- 





Lagrange's algorithm 
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We now combine this with the results of the foregoing section in order to obtain a 
coherent description of the computational procedure to determine the li, Xij. 

When oi, 02, 03 are given, one first has to set up the basis (ei, e2) of L given in Lemma 
4(a). This requires 0(log(maxi Oi)) steps for the euclidean algorithm to find mi, m2. 
Using this basis as input one then runs Lagrange's algorithm and obtains a reduced 
basis (u, v)< of L. As shown in [3], Theorem 3.0.3, this is accomplished in at most 
0(log(maxi Ci)) loops of the algorithm. Replacing u by — u, if necessary, makes Ui > 0, 
Uj < (j 7^ i) for some i. Now put A := Vi/ui and define v_, v_|_ as given in Theorem 
7. The quantities li, xij can then be read off from the coordinates of the three vectors 
u, v_, v_|_. Note however that the bit complexity of the whole computation is quadratic 
insofar as the bit-length of the input data a^ not only affects the number of required 
steps/loops but also the cost of each step. 

Example 2. Consider the random numerical example oi := 4327, 02 := 6716, 03 := 9237. 
In 8 steps the euclidean algorithm finds 2055ai — 1324a2 = 1, and after 6 loops of 
Lagrange's algorithm we arrive at a reduced basis of L given by u = (—53,-47,59), 
V = (-130, 59, 18). Theorem 7 then tells us that 

fi = (130,-59,-18), £2 = (-77,106,-41), £3 = (-53, -47,59) 

is a solution basis for the given data, and by Theorem 3 we have 

5(a) = 920 947 , iV(a) = 493 045 . 
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